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Abstract 

Out of self-consistent semi-classical calculations performed within the so-called Ex- 
tended Thomas- Fermi approach for 212 nuclei at all even angular momentum values 
/ ranging between and 80 h and using the Skyrme SkM* effective force, the I- 
dependence of associated liquid drop model parameters has been studied. The latter 
have been obtained trough separate fits of the calculated values of the strong inter- 
action as well as direct and exchange Coulomb energies. The theoretical data basis 
so obtained, has allowed to make a rough quantitative assessment of the variation 
with I of the usual volume and surface energy parameters up to spin of ~ 30-40 H. As 
a result of the combined variation of the surface and Coulomb energies, it has been 
shown that this /-dependence results in a significant enhancement of the fission 
stability of very heavy nuclei, balancing thus partially the well-known instability 
due to centrifugal forces. 
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1 Introduction 



The availability of heavy ion beams at sufficiently high energy has made 
it possible to transfer very high angular momenta to nuclei, particularly in 
fusion-evaporation channels. Upon increasing the angular momentum one has 
opened up the exploration of a wide range of very large nuclear deformations, 
as in superdeformed states which were previously reached only in the fission 
process of ground (or weakly excited) states of heavy nuclei. As a matter of 
fact, one is currently able to reach the critical angular velocity regimes where 
the nuclear system experiences a centrifugal disassembly on all the accessi- 
ble part of the nuclear chart. This is why a long time ago already, various 
authors have studied the equilibrium shapes and stability properties of ro- 
tating nuclei within a liquid drop model approach, drawing a close analogy 
between the properties of such a self-bound mesoscopic system with those of 
self-gravitating and rotating celestial objects (see Ref. [1] and Refs. quoted 
therein). Later, a more microscopic description of such nuclear systems has 
been provided upon minimizing the so-called "Routhian" variational quantity 
in a macroscopic-microscopic approach (from the pioneering works of Refs. 
[2,3] to more systematic approaches as e.g. in Ref. [4]) thus relying again on 
the liquid drop model assumption when using in practice the Strutinsky energy 
expansion [5] in a slightly, and easily, transposed version. More recently, fully 
self-consistent Routhian approaches (after preliminary attempts making use 
of an inert core approximation [6,7]), have been made available either with the 
Skyrme [8-11] or the Gogny [12,13] parameterization of the nucleon-nucleon 
effective interaction and have been shown to be reasonably realistic. 

While the latter type of calculations does not imply any liquid drop model 
assumption to study the behavior of rotating nuclei, they are not really easy to 
handle in systematic calculations and certainly not suited, due to shell effects, 
to provide a transparent account of the bulk properties of nuclei at finite 
angular momenta. Yet, one would like to take stock on the impressive ability 
of such effective force parameterizations (Skyrme SkM* [14] or SLy4 [15] and 
Gogny D1S [16] interactions) to yield a good reproduction of a wide array of 
static and (low excitation energy) dynamical properties. As well known, the 
theoretical framework to search for the bulk nuclear properties embedded in 
such purely microscopic approaches is a semi-classical approach a la Wigner 
[17], i.e. performing a truncated expansion of the solution (namely of the one- 
body reduced density matrix in the Hartree-Fock case). More specifically, we 
will use here the Extended Thomas-Fermi approach proposed in Ref. [18] and 
exhaustively discussed in Ref. [19]. 

In our Routhian approach we make use of a Skyrme type effective force pa- 
rameterization. We explicitly express the functional dependence of some local 
density functions (as e.g. current and spin-vector densities) providing the rel- 
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evant informations on the non-local (in r) part of the time-odd component of 
the one-body density matrix, in terms of the spin-scalar density p(r). Some 
steps in that direction have been made by Grammaticos and Voros both with 
[20] and without [21] spin degrees of freedom. However the full analytical 
solution of this problem has only been provided recently [22] and used to per- 
form self-consistent semi-classical calculations [23] . The present description of 
bulk nuclear properties at finite angular momentum is clearly an extension of 
these previous studies. It fully develops some ideas presented earlier [23] in a 
somewhat exploratory fashion. 

Indeed, we want to assess the relevance of the concept of a liquid drop for the 
description of rotating nuclei at equilibrium. The rationale for doing it stems 
from the anisotropic character of the rotational constraint imposed on the 
considered piece of nuclear matter which seems at first sight to be at variance 
with the isotropic character attached to a saturated liquid drop. Furthermore 
nuclei being leptodermous, as nicely illustrated by Myers and Swiatecki [24], 
one may wonder whether or not the surface tension, and hence the skin prop- 
erties, might not be affected by the rotational character of the nuclear fluid. 
The above may be summarized in a practical form by recalling first that in 
Ref. [1], for instance, the nuclear energy is written, with usual notation, as 

E(P,I) = E(P,I = 0) + ^, (1) 

where (3 stands for a set of deformation parameters. One takes now for the 
energy E((3,I = 0) a standard version of the liquid drop model and the rigid 
body moment of inertia associated with the given shape for the parameter 
^ 2 K Therefore one assumes that all the /-dependence is only contained in the 
kinetic collective energy of a rigid rotation. 

Let us translate that into the Hartree-Fock language. One may then split the 
total energy in the laboratory frame into two pieces: 

E = E + + E_, (2) 

the first involving only the part p + of the density matrix which is even with 
respect to the time- reversal symmetry and the second, E_, issuing from the 
time-odd part p_ of the density matrix (despite the two-body character of the 
interaction part of the Hamiltonian H, there are no cross-terms in Eq. (2) 
due to the time reversal invariance of H). For the standard Skyrme energy 
density functionals, the former part involves local densities (for both charge 
states) like the spin-scalar diagonal densities, the kinetic energy densities and 
the so-called spin-orbit densities (see e.g. Ref. [25]), while the latter involves 
the current densities and the spin- vector diagonal densities (see e.g. Ref. [9]). 
In the adiabatic limit, such a partition into two energies does correspond to 
the partition assumed in Eq. (1). In non-adiabatic regimes, the dependence of 
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the two parts p + and p_ of p in terms of the collective velocity (in the present 
case the angular velocity a;) is more complicated since p + may contain, for 
instance, a quadratic dependence on uj. 

One is thus lead, in practice to the two following questions: 

- Is it correct that in a self-consistent semi-classical approach E_ may be 
written as in Eq. (1) with a rigid body moment of inertia? 

- Is it correct to assume that E + may be well approximated by a liquid drop 
energy, and if so are the parameters of this liquid drop independent on the 
angular momentum? 

The answer to the first question is essentially affirmative as it has been demon- 
strated in Ref. [22]. There, it has been shown that the semi-classical moment 
of inertia S may be split into two parts for which compact mathematical ex- 
pressions have been given in the Skyrme force case. The first one is the zero 
order term of the semi-classical expansion (the so-called Thomas- Fermi term) 
and is as expected (and well known before the work of Ref. [22]) the rigid 
body moment of inertia. The remaining part corresponds to the second order 
term of this expansion which can be decomposed into two contributions, one 
coming from the orbital, the other from the spin degrees of freedom. They 
may be illustrated upon using a well established analogy between an angular 
velocity and a magnetic field B (analogy between, for instance, the Coriolis 
pseudo-force and the Lorentz magnetic force, first advocated in the context 
of a rotating many fermion system by Dabrowski [26]). For realistic effective 
forces (in Ref. [22], we have used the SkM* force which is known to give, at 
least in the /^-stability valley a very satisfactory account of bulk nuclear prop- 
erties) one has found a negative contribution of the orbital degrees of freedom 
to corresponding to an analog of a Landau diamagnetism, while the spin 
degrees of freedom contribution being positive, corresponds to a Pauli para- 
magnetic alignment. Both contributions are rather weak with respect to the 
Thomas-Fermi contribution, and moreover, as we have seen, they cancel each 
other to a large extent. Specifically, it has been found in Ref. [22] that one can 
write with obvious notation for a nucleus of mass A 



with the following values for the coefficients of the second order contributions 
(again for the SkM* force): 



leaving thus only a small correction to the rigid body ansatz, especially for 
heavy nuclei. 

The problem of the kinetic collective energy being settled, the second question 



(3) 



rj t m -1.94, r] s m 2.63, 



(4) 
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raised above is still open. First, is there any /-dependence in the non-kinetic 
(E + ) part of the energy? The answer to this is clearly yes. In Ref. [23], we have 
found that the self-consistent semi-classical energy E + , for a given nucleus and 
using the SkM* force, was indeed very much varying. Therefore, if at all valid, 
the liquid drop model to be used should clearly incorporate coefficients varying 
with the rotational velocity. An important question needs to be clarified at 
this point. A priori, there are no clear reasons why the liquid drop parameters 
should vary with / or uj . Of course, the former being an observable which is 
furthermore a conserved quantity, seems a better candidate. But it is not clear 
to us why this quantal feature should reflect itself in this way for the rotational 
dependence of a rather crude modelization of the energy. However in Ref. [23], 
it has been shown that the variation of the energies E + were of similar relative 
magnitude between light and heavy nuclei, only when compared for the same 
value of / (and not for the same value of uo which scales as A~ 5 / 3 for a given 
value of I). We are thus lead to conclude phenomenologically that, if valid, 
the rotating liquid drop model should be /-dependent and not w-dependent. 
(Note that even though, as stated above, the results of Ref. [23] may lead to 
such a conclusion, the liquid drop model parameter fit in terms of the angular 
velocity performed there is somehow inconsistent). 

The remaining issue is then the validity of the liquid drop model at finite 
angular momenta. For that purpose, we need to perform for each value of the 
angular momentum which we decide to study, high accuracy self-consistent 
semi-classical calculations for a sufficient number of nuclei covering reasonably 
well the chart of nuclides so that we should be able to make a significant 
liquid drop parameter fit. This was not at all the case of the calculations of 
Ref. [23] (too few nuclei, 33, and an accuracy in the variational semi-classical 
solutions somewhat insufficient for the present purpose). It is the aim of the 
present paper therefore to achieve in a numerically satisfactory fashion the 
preliminary calculations sketched in Ref. [23]. 

Two remarks are worthwhile at this point. First, it is clear that we do not 
include pairing correlations. This could be done, as we are currently attempt- 
ing, within the BCS or the Hartree-Fock-Bogoliubov framework in the semi- 
classical approximation (for the Thomas-Fermi and Extended Thomas-Fermi 
non- rotating case see Refs. [27] and [28] respectively). The second remark deals 
with the bulk deformation properties. All our self-consistent semi-classical cal- 
culations will be performed at sphericity. As a matter of fact when assessing 
the validity of a liquid-drop model approximation, one intends to do two dif- 
ferent things. The first is to check the ability of a liquid drop parameterization 
to reproduce with a good accuracy the semi-classical energies at a given defor- 
mation. For that purpose choosing spherical solutions is an obvious choice for 
the sake of computational convenience. The second is to check whether or not 
the semi-classical energies are varying with deformation like a sum of various 
geometrical quantities associated with the deformation (surface, Coulomb en- 
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ergy or more). This will not be attempted here. Yet, we think obvious that the 
present approach casts some light on nuclear deformation properties at finite 
angular momentum (through the fissility parameter for instance). Restricting 
ourselves to spherical shapes, we may therefore ignore the difficulties described 
by Krappe, Nix and Sierk [29], and studied in the rotational context in Refs. 
[30,31], when dealing with nuclear shapes having some concavity (like fission 
saddle points of some nuclei) or not simply connected (like post scission or 
fusion configurations, for instance). 

The paper will be organized as follows. Section 2 will detail the theoretical 
framework in which we have studied the bulk rotational properties of the nu- 
clei. In Section 3 we will present the calculational approach insisting on the 
choice of the sample of nuclei, the accuracy of the self-consistent calculations 
and the fitting procedure. The results will be displayed and discussed in Sec- 
tion 4. Finally, Section 5 will contain some conclusions and perspectives of the 
present work. 



2 The theoretical framework 

To describe in a semi-quantal way the wave function of a nucleus rotating with 
respect to the laboratory frame with an angular velocity ou, one minimizes the 
Routhian R defined in terms of the microscopic Hamiltonian H as 

R = H-uI. (5) 

In the Hartree-Fock approximation, the single particle states defining the solu- 
tion is thus reduced to an eigenvalue problem involving the one-body reduction 
of the Routhian of Eq. (5). The corresponding semi-classical one-body density 
matrix may thus be obtained according to the usual Extended Thomas- Fermi 
method upon replacing the Hartree-Fock one-body Hamiltonian by the cor- 
responding one-body Routhian. For a Skyrme effective force, the Routhian 
expectation value is obtained as an integral over space of a Routhian density 
which is expressed (see e.g. Ref. [9]) in terms of various densities (spin-scalar 
and spin- vector densities, kinetic energy densities, spin-orbit and current den- 
sities). In the semi-classical Extended Thomas- Fermi approach all densities 
may be expressed as functionals of the spin-scalar densities [22]. Therefore 
the Routhian density is itself a functional of the spin-scalar densities. The 
variational solution for a given angular velocity may be just obtained then by 
minimizing the integral of this Routhian functional. 

To do so, choosing the Oz the rotation axis, we have used a parameter- 

ization of the density functions p q (r) (here called p^ ew for reasons explained 
below), where q stands for the charge state, expressed in cylindric coordinates 
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(r, z, (p) as 



P ; icw (r, z) = p 0q [l - a q exp(-/3 9 2 r 2 S f 
1 + exp 



r 2 Sq Z + Z 2 S g — R q 



-7q 



(6) 



Whereas we want to describe spherical solutions, they appear a priori to have 
only a cylindrical symmetry (see however the discussion below). The densities 
depend on 14 parameters (7 for each charge state). However the space integrals 
of the densities p g (r) should be equal to N q , the particle number of charge q, 
leaving thus 12 independent parameters. For a q = 0, and s q — 1, one is left 
with surface asymmetric (for 7 g ^ 1) Fermi distributions as widely used in 
Ref. [19] for instance 



Pq{r) = POq 



1 + exp 



' y/x* + y 2 + z* - R q 



(7) 



One notices in Eq. (6) that in front of the modified Fermi term we have en- 
closed a factor similar to what has been considered in Ref. [19] with a very 
significant difference, namely its anisotropic character. This has been intro- 
duced to let the solution adjust to the centrifugal forces which are themselves 
not isotropic. For the same reasons advocated in the slightly different context 
of Ref. [19], we have reduced the variational space by imposing that f3 q should 
have a fixed value within a range where this modification would not alter the 
desired leptodermous and smooth character of the semi-classical densities. In 
practice we have used: 



(3 2 q = l/{c q R q )\ c q = 0.5, 
Rq being the radius parameter introduced in Eq. (6) 



(8) 



Let us discuss now the scaling parameter s q appearing in Eq. (6). Introducing 
an anisotropy as above described, our solution lacks its spherically symmet- 
rical character. To be able to compare solutions with various values of the 
parameters a q we must constrain them to keep, in some way, a vanishing 
deformation. This is achieved by imposing that 



(9) 



where {r 2 ) q and (z 2 ) q stand for the expectation values of r 2 and z 2 for the 
distribution of nucleons of charge q. Taking into account the particle number 
conservation, one gets the following explicit expression for the parameter s q : 



roo re 

5~ 2 = 2 / r dr 
H Jo J— 



z 2 p° ld dz 



dr 



PT dz 



(10) 
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with 

'Vr s Tz s -R Q ^' l ~~ ! " 



P° q ( r > z ) = Pog 1 - "<? exp(-/3 r ) 



1 + exp 



Oq 



:n) 



To obtain the expression for s 2 one may introduce stretched coordinates r = 
s^ 3 r and z = s~ 2 ^ 3 z in terms of which one has p q ew (f, z) = pf d (r, z). 

The usual iterative procedure now begins with the determination of the scaling 
factor s q . Then we find the parameter p 0q from the normalization condition 



4tt / r dr / (7?™{r,z) d ^ ( 12 ) 
jo Jo 

Hereafter JV p = Z, N n = N. In all the integrals we take advantage of the 
reflexion symmetry p q ew (r, —z) = p q cw (r, z) . 

Once the variational parameters have been determined, one obtains the total 
energy in the laboratory frame according to the decomposition of Eq. (2) 
where E + is obtained by using the time-even energy Skyrme functional of Ref. 
[25]. One gets namely (see Ref. [22]). 

E[p] = J M[P\ d 3 r = / 3WM d 3 r + \& 2 \p)u 2 , (13) 



with 



^statb] = ^ £ /^(P) + C ntml (P) + ^(P). (14) 



2m q 



where f q is the usual Skyrme effective mass factor and r q (p) the ETF expres- 
sion for the static kinetic energy density up to /i 4 -terms (see e.g. Ref. [19]). The 
expression for V q central (p) is given with the notation of Ref. [9] for the Skyrme 
parameters (Bi, a) where as a general rule all densities without a charge state 
subscript q are total densities (i.e. for instance p = J2 q Pq)- 

KT tml = B lP 2 + B 2 J2p 2 q + B 5P V 2 p + B 6 £ Pq V 2 Pq 
i i 

+ (B 7 p 2 + B s J2p 2 q )p« (15) 



q 



The ETF expression for the spin-orbit contribution up to h -terms V q °' to 
^stat(p) is given e.g. in Eq. (A5) of the review paper by Brack et al. [19]. 
Finally, the ETF expression of Q^ 2 ) up to fr 2 -terms has been derived in Ref. 
[22] (see e.g. Eq. (63) there). 

To obtain solutions corresponding to precise values of the angular momentum 
J, we adjust iteratively the value of u accordingly. We have made calculations 
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Table 1 

Isotopes used in the fit. For each elements we have considered all even-A isotopes 
from Ai to Af (including the boundaries) corresponding to a number of isotopes 
N is . 



Element 


Z 


Ai 




N is 


Element 


z 


M 




N is 


Ca 


20 


34 


78 


23 


Yb 


70 


166 


180 


8 


Ni 


28 


56 




1 


Os 


76 


184 


198 


8 


Ge 


32 


58 


118 


31 


Pb 


82 


182 


242 


31 


Zr 


40 


86 


100 


8 


Th 


90 


226 


240 


8 


Sn 


50 


94 


168 


38 


Pu 


94 


236 


252 


9 


Ce 


58 


134 


148 


8 


Cf 


98 


252 




1 


Sm 


62 


120 


194 


38 













for a sample of 212 even-even nuclei, covering all the chart of nuclides whose 
mass is higher than 34 within reasonable estimates of the two drip lines and 
including five long isotopic series (Ca, Ge, Sn, Sm, Pb), as detailed in Table 
1. Obviously for semi-classical solutions whose static properties are smooth 
functions of N and Z, the particular choice of nuclei is irrelevant provided that 
the sampling is dense enough in the region of interest for the properties under 
study. For instance to get reasonable estimates of proton-neutron asymmetry 
properties it is essential to include in the sample isotopes up to the drip lines. 
We will discuss in the next Section the quality of our sample in this respect. 

Now, for each value of J, we have made a liquid drop model fit of the 212 
nuclear energies E nuc , defined in terms of the energy E + as 

-Enuc = E + — E cou \, (16) 

where E cou \ is the total Coulomb energy. The latter includes the direct term 
(.Ecd) and the so-called "Slater approximation" (namely a local density ap- 
proximation [32] of the Bethe-Bacher [33] nuclear matter estimate) for the 
exchange term (Ece) whose accuracy has been checked earlier [34]. Both 
Coulomb energy terms have been fitted separately (i.e. independently of each 
other and of the nuclear energy). We leave for another publication [35] an in- 
depth discussion of the liquid drop fit of the Coulomb energies, in particular its 
isospin-dependence which is generally overlooked (see a preliminary account 
of such a work in the review of Ref. [36]). Namely, the liquid drop expression 
for E nuc which we have considered here, is given for the nucleus labeled by j 
(j — 1, • • • Na, where Na is the number of nuclei entering the fit) as 
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EUl) = a v (I) 
+a s {I) 



l-k v (I)(l-2Z j /A j f 
l-k s (I)(l-2Z j /A j 



A, 



1 2/3 



(17) 



The Coulomb part is the sum of a direct (Ecb) and exchange (Eqe) part 
whose absolute value is known to be one order of magnitude smaller than Ecn 
(see e.g. Ref. [34]). For the direct term we have made two fits one which is 
rather usual (sharp density term plus a leptodermous surface correction as in 
Ref. [37]) 



EUn = c(I) ZVAf-^1) Z*/A„ 



(18) 



and another extended one discussed in a forthcoming paper [35] which takes 
into account both the surface diffuseness and (N — Z) isotopic effects 



EUl)=c(I) 

-7(/) 



1 + ai(/)(l - 2Z J /A j )Aj 1/s + a 2 {I)Aj 1/3 
l + P(I)(l-2Z j /A j )] Z]/Aj. 



Z)jAf 



(19) 



For the exchange term (-Ece) we consider the usual exchange term deduced 
from an infinite nuclear matter non-diagonal density matrix [33] 



EUI) = c'(I) zflAf. 



(20) 



3 Some computational details 



We use the SkM* force with the original set of parameters previously de- 
termined by fitting HF results to experimental data. The fit of the effective 
interaction has been made correcting approximately the HF energy by sub- 
tracting only the one-body part of the mean energy of the center of mass (cm.) 
motion, (see Sect. 2.3 of [38]). To be consistent, in our ETF calculations 
we will also only take into account that part of the cm. energy correction. The 
introduction of this correction leads to the replacement of the nucleon mass 
m by mA/{A — 1) in the total kinetic energy and so in our non-local case we 
have to replace Y, q f q T q (p) by ^2 q f q (l — \/A)r q {p) in the energy density 
[see Eq. (14)]. 

The calculation of the Coulomb energy is the most time consuming part of the 
problem due to its range. We take into account both the direct and exchange 
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parts of this energy: 

£coui = \ J d 3 n p P (r 1 )Vc(r 1 ) + E°£ uh (21) 

with 

Vb(rx) = e / d 3 r 2 ^£L. (22) 

The Coulomb exchange part is treated as already mentioned in the so-called 
"Slater approximation" [32-34]: 

E^a = ~(3M 1/a Jd i rp^(r). (23) 

The integrand in the direct part of the Coulomb energy has a logarithmic 
singularity at the point r*i = r 2 . Therefore such a formula is not suitable 
for numerical integration. A way to bypass the latter difficulty is to use the 
relation 

V^|r 1 -r 2 | = 1 , (24) 

|ri-r 2 | 

and to carry out two integrations by parts to get: 

Vein) = e -J d 3 r 2 |n - r 2 \ V 2 r2 p p (r 2 ). (25) 

With our parameterization (6) of the density it is most convenient to treat the 
problem in cylindrical coordinates choosing the axis of symmetry to be the 
Oz axis. Integrating over the azimuthal angle <p 2 the expression for Vc(j\,zi) 
reduces to [39] 

rOQ roo I 

Vein, zi) = 2e r 2 dr 2 / d^ 2 J (z 1 - z 2 f + [r x + r 2 ) 2 

JO J-oo v 

E(k)V 2 p p (r 2 ,z 2 ), (26) 

with 

k 2 = 4nr 2 [(^i - z 2 f + (n + r 2 ) 2 ] 1 (27) 
and where E(k) denotes the complete elliptic integral of the second kind: 

E(k) = Jl -ksin 2 u du. (28) 
Jo 

We have therefore for the direct part of the Coulomb energy: 

^Coui = ™ / n dn / dzi p p (rx, ^^(n, Zi). (29) 

JO J-oo 
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In practice, the infinite limit of integration is replaced by R = 25 fm, and we 

get 




f(r 1 ,z 1 ,r 2 ,z 2 ) + f(r 1 ,z 1 ,r 2 ,-z 2 ) . 



(30) 



This 4-fold integral is evaluated numerically using the routine dOIfcf of the 
NAG library. The same code has been used for the computation of the integral 
expressing the number conservation condition. In the simpler case of a q = 
and s q — 1, i.e. for the parameterization (7) of the densities, a Gauss-Legendre 
quadrature formula has been used in both integrals (where the integration 
interval has been divided in 4 equal pieces in which 16 integration points have 
been considered). 

In calculating the rest of the energy we use a trapezoidal rule with a step 
of 0.05 fm. For that part, the integration both in r and z coordinates goes 
up to R = 24 fm. Different tests on the precision of the ETF calculations as 
connected to the used procedures for numerical evaluation of the integrals have 
been made, of which we will merely give here some examples. In particular, to 
illustrate the quality of our estimation of the direct Coulomb energy as given 
by equation (29), for the case of the heavy 208 Pb nucleus, we may mention 
that upon computing it with the same mesh size for R = 25 fm and R = 50 
fm we obtain the same value up to seven digits (namely up to better than 
100 eV). The relative accuracy for the calculation of E nuc is of the order of 
10~ 6 . For instance upon increasing the mesh size by 25 percent and using the 
standard above described cut-off, the nuclear energy of 208 Pb varies relatively 
by less than 10~ 7 . Conversely, for the mesh size in use in our calculations when 
increasing the cut-off value of R = 24 fm by 25 percent, the corresponding 
variation for the same 208 Pb nucleus is again of the order of 10~ 7 . For the 
lightest nuclei under consideration here it is one order of magnitude larger 
(but still better than 1 keV in absolute terms). 

The optimization procedure of the density has also been checked. The termi- 
nation criterion is delicate in any multidimensional minimization routine and 
in the case of the density of Eq. (6) we have 8 independent parameters to be 
optimized. In practice we use a downhill simplex method (see e.g. [40]) and 
the numerical code stops when the decrease in the function value in the termi- 
nating step is fractionally smaller than some tolerance. In our calculation this 
tolerance was set to 10~ 7 . Such criteria might be fouled by a single anoma- 
lous step. In our case however, we may use the fact that the calculated ETF 
energies are semi-classical and thus, all the parameters are expected to vary 
slowly when changing the number of particles, so that when the minimization 
is thoroughly performed for one nucleus, we know that the deviation of the re- 
sults for a neighbouring nucleus should be rather small. The same is expected 
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when going from a given value of the angular momentum to a neighbouring 
one (i.e. here, differing by 2 h). Therefore to make sure that we have reached 
a good local minimum, the calculations have been run using as initial values 
of the parameters at a given angular momentum, the optimized parameters 
as obtained at the preceding value of the angular momentum (i.e. in our case 
lower by 2 K). To illustrate this smooth behaviour, we display on Figs. 1-3, 
the angular momentum dependence of the density parameters for four nuclei 
90 Zr, 150 Sm, 208 Pb and 240 Pu. For the three last nuclei, one gets gently varying 
density parameters up to spin values of the order of 60 h. The overall quality 
of these parameters is slightly worse for 90 Zr. 

Let us go back to the tolerance value adopted for the termination of the 
density parameters optimization. Such a level of relative accuracy has been 
found necessary in order to ensure the stability of the optimization process. 
However, it has been above noted that for light nuclei the computation of the 
integral with respect to the cut-off radius was one order of magnitude larger. 
In doing the optimization we must therefore assume that the latter numerical 
error in evaluating the energy will not affect significantly the localization of 
the optimal parameters. This has been shown to be true in a limited number 
of test cases at zero spin. 

Let us now briefly describe the fitting procedure in use here for the energy. 
Mathematically, one minimizes for each value of the angular momentum Ih 
the weighted least-squares function 

Na r i2 

xV) = Y,^[EfV)-E 3 Mi))] , (31) 

with £^ alc (7) being the nuclear part of the ETF calculated energies or the 
Coulomb direct /exchange energies, Ej(a,k(I)) the corresponding expression 
with p parameters a,k(I) to which we want to fit the calculated energies and 
the weights Wj have all been set to 1. As before, Na is the number of nuclei 
used in the fit. 

For this fit we used the code afxy [41] which is an interactive numerical re- 
searcher of nonlinear equations and nonlinear least square fits. This algorithm 
is based on a composition of autoregularized Gauss- Newton methods with im- 
proved versions of the Levenberg-Marquardt method and the Tikhonov-Glasco 
descent-on-parameter method [42] . However, our fitting problem is rather sim- 
ple and it appears that a Gauss-Newton iteration process is sufficient. The 
resulting root mean square deviation between calculated and fitted energies, 
is defined as 



AE 



N A 



Y,[EfV) - Ejiafil))] , (32) 
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Table 2 

Fitted LDM-parameters and estimated errors (absolute values) for different spin 
values. 



i/h 





20 


40 


60 




fit error 


fit error 


fit error 


fit error 


a v 
k v 


-15.453 0.007 
1.748 0.009 


-15.444 0.007 
1.743 0.008 


-15.436 0.007 
1.743 0.008 


-15.472 0.013 
1.775 0.016 


a s 
k s 


16.961 0.037 
2.020 0.041 


16.908 0.036 
1.999 0.040 


16.848 0.037 
1.996 0.041 


17.016 0.069 
2.144 0.075 



where a£ p (/) are the optimized parameters (i.e. obtained upon minimizing the 
X 2 (/)). Standard deviations of the parameters have been obtained by usual 
numerical methods. The corresponding errors are listed in Table 2, together 
with the corresponding fitted values, for the four parameters (a v , k v ,a s ,k s ) of 
E imc at four different values of the angular momentum. As a result it appears 
that the significance of the fitted values is much lower above I — 40 h, while 
it remains qualitatively comparable below this value (we refer to the next 
Section for a further discussion of the quality of the fit in terms of the root 
mean square energy deviation). 



4 Results and discussion 

At first, we discuss the usefulness of the degree of freedom introduced in the 
density to allow for a spatial anisotropy with respect to the angular velocity 
[see Eq. (6)]. As a matter of fact, upon performing two series of calculations 
(constraining the a q parameters to vanishing values or letting them free in 
the minimization procedure) we found, up to rather high values of uo (up to 
huj = 1.25 MeV for a spherically constrained 208 Pb nucleus corresponding to 
/ ~ 130 h), that the gain in energy due to this degree of freedom was small 
(~ 110 keV) and constant (up to a couple of keV). We therefore chose to 
ignore this effect in what follows (i.e. set a q — 0). 

We will first consider the strong interaction part of the nuclear energy. Let us 
assess the quality of our fit in terms of the root mean square energy deviation 
AE. For that we will study AE as a function of a lower bound value A m - m 
of the nucleon number of the retained nuclei (A > A m i n ). The corresponding 
number of nucleons are reported in Table 3. 

The resulting energy deviations are plotted on figure 4. As a result, it appears 
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Table 3 

Number of nuclei considered in the fit for a given lower bound A m [ n of the total 
nucleon number Annuel- 



A ■ 


34 


40 


50 


60 


70 


80 


90 


100 


110 


Annuel 


212 


209 


204 


197 


181 


177 


170 


157 


146 



that for angular momenta which are not too high (lower than 30-40 H) and for 
intermediate values of A min (in the 60-80 range) the quality of the fit seems 
rather stable. For large values of A min the size of the data basis is probably 
too limited. On the other hand light nuclei at finite angular momenta / seem 
to be not adequately represented by a liquid drop formula (as seen on Fig. 4). 
This is more and more so upon increasing /. 

As another test of the consistency of our data basis and of its physical limita- 
tions, we have studied the convergence of the fitted LDM parameters a v and 
a s (Fig. 5) as well as k v and k s (Fig. 6) in terms of the number of considered 
nuclei Np varying as described below. We define an approximate valley of sta- 
bility as a function N(Z). The minimal distance of each considered nuclei to 
this valley is then computed. Increasing Np corresponds to including into the 
fit nuclei which are more and more distant from the stability curve. 

The results displayed on Fig. 5 clearly demonstrate that the data basis size 
is just barely adequate for rough estimates of the dependence of the a v and 
a s coefficients as functions of I at least for moderate spin values (/ < 30-40 
h). As seen on Fig. 6, at the present stage, the number of retained nuclei does 
not allow to draw safe quantitative conclusions on the /-dependence of the 
asymmetry coefficients k v and k s . 

The results obtained on Fig. 5 are consistent with a quadratic dependence of 
the a v and a s parameters as 

a v (I) = a v (0) (l + c v I 2 ) (33) 

and 

a s (I) = a s (0) (l + c s I 2 ) . (34) 

The results of the fit of c v and c s , upon using all a v (I) and a s (I) calculated 
for I going from to A0h by 2 ft-steps, are plotted on Fig. 7 as functions 
of the number of considered nuclei Np (defined as above described). These 
parameters are not converged, yet up to about 10 percent one may assign to 
them the following "asymptotic" values 

c v « 2.0 10~ 5 and c s « 1.1 10~ 4 . (35) 



15 



Let us now discuss the variation with angular momentum of the Coulomb 
energies. We present in Fig. 8 the evolution of the root mean square radius of 
the proton density as a function of the angular momentum for the four already 
considered nuclei. Due to the effect of centrifugal forces, it is expected that the 
radius should increase with an increasing angular momentum which is indeed 
the case here. As a consequence of their inverse proton radius dependence, the 
Coulomb energies (both direct and exchange terms) plotted in Fig. 9 exhibit 
decreasing patterns in absolute values as functions of I. 

We have then performed a fit of the direct Coulomb energy E CD using the 
two-parameter formula of Eq. (18). Since the Coulomb energy is a decreasing 
function of /, one may have expected also to observe a decreasing value of 
the parameter c corresponding to the leading contribution to Ect>- As shown 
on Fig. 10 this does not turn out to be the case. This incoherence is, in our 
opinion, due to the inadequacy of Eq. (18) to describe consistently the isospin 
dependence of the Coulomb energy. As described in Ref. [35], the 2-parameters 
formula of Eq. (18) provides a very poor fit of Eqb- Indeed it yields at zero 
spin for the same data basis (i.e. including nuclei up to the drip-lines) a RMS 
energy deviation of some MeV, which is close or above the maximum of the 
expected spin-dependence of E CD . This is why we have performed the fit with 
the 5-parameters formula of Eq. (19). With this new parameterization, we 
obtain a parameter c that now decreases with increasing I (as shown on Fig. 

11 together with the variation of the other fitted parameters). 

As for the fit of the exchange term Ece, we obtain as demonstrated on Fig. 

12 a d parameter which decreases as a function of the angular momentum. 
This behavior is consistent with the above described variation of the proton 
radius. 

The variations with I of the parameters of both the strong interaction and 
the Coulomb part of the liquid drop model energy entail some important con- 
sequences on the nuclear deformability and hence on collective excitations as 
well as on the fission stability. To illustrate the latter point we have evaluated 
the fissility parameter x defined as usual by 

F° 

x = =^ (36) 

(where Eq ou1 and E® stand for the Coulomb and surface energies at zero defor- 
mation) of the 236 U nucleus for spin values varying from to 40 h. In doing so, 
we took the surface energy as well as Coulomb direct (5 parameters fit) and 
exchange energies resulting from our fit. It is to be noted that even though the 
variation with / of the parameters k s is not very well known quantitatively, the 
corrective character of the asymmetry term with respect to the total surface 
energy (~ 10 percent) allows us to draw anyway significant conclusions on the 
variation of x with /. As shown on Fig. 13, the fissility parameter decreases 
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by ~ 2.5 percent when increasing the angular momentum from to 40 H. 
According to the liquid drop estimate for the fission barrier height Bf of Ref. 
[43] 

B f ~ 0.83 E° s (1 - x) 3 (for x > 2/3) (37) 

One sees that such a variation of x corresponds to an increase of Bf by approx- 
imately 30 percent for the considered 236 U nucleus. This result, far from being 
insignificant, should be put into the perspective of the rotating liquid drop 
approach of Ref. [1]. Upon increasing the angular momentum, the rotational 
energy induces as well known an increased instability against fission. The vari- 
ation of the fissility parameter with the spin which we have demonstrated here 
reduces the effective amount of the resulting fission instability. 



5 Summary and conclusions 

In this work, we have attempted to estimate the angular momentum depen- 
dence of the liquid drop model parameters and hence the validity of so-called 
rotating liquid drop model approaches making use, as well known, of a spin- 
independent liquid drop energy beyond the usual rotational energy term. To 
do so, we have first performed for 212 different nuclei at 41 different angular 
momentum values, state of the art self-consistent semi-classical calculations. 
Their results correspond to the h — > solutions (beyond the trivial Thomas- 
Fermi approach) of a time-reversal breaking Routhian problem using the full 
Skyrme SkM* effective force. A particular attention has been devoted to the 
numerical stability of the corresponding variational solutions. The obtained 
strong interaction and Coulomb energies have been separately fitted into liq- 
uid drop model expressions at each angular momentum value. As a result, we 
have shown that our theoretical data basis was just barely capable of providing 
rough estimates of the angular momentum variations of the standard volume 
(a v ) and surface (a s ) parameters while providing only qualitative insights into 
the behavior of the asymmetry parameters (k v and k s ). The combined variation 
of the surface and Coulomb energies reflects itself into a sizeable variation of 
the fissility parameter which yields, as a consequence, an important enhance- 
ment of the fission stability for very heavy nuclei, thus partially balancing the 
well-known instability generated by the centrifugal forces. 

Clearly, one possible continuation of the present work should consist in con- 
siderably enlarging the theoretical data basis in such a way as to extend the 
level of confidence in the numerical estimates of the resulting c v and c s pa- 
rameters assessing the variations according to the angular momentum of a v 
and a s . Another goal should be to have also access to the /-dependence of k v 
and k s . The interest of such a numerical study could be balanced by the fact 
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that as soon as one is only interested in the behavior of a few nuclei, one could 
as well perform "exact" self-consistent semi-classical calculations. However it 
could be argued that such calculations would necessitate, when performed for 
very deformed nuclei, or pieces of nuclear matter exhibiting crevices or not 
simply connected, etc., highly non trivial generalizations of the model densi- 
ties in use for the variational problem. At any rate, it remains that our present 
study clearly points out the necessity to go beyond the use of a zero angu- 
lar momentum parameter set when dealing with a rotating liquid drop model 
approach. 

Some other lines of further work could be mentioned. One would consist in 
evaluating carefully and systematically the combined influence of the rota- 
tional energy term and of the spin-dependent liquid drop energy term on 
fission properties. The other set of fascinating questions stems from the neg- 
ative result obtained here when trying to identify an anisotropy effect in the 
spatial matter distribution. It is indeed not clear to us that such a negative re- 
sult should survive whenever one would deal with strongly deformed solutions. 
Moreover, the semi-classical solutions in use here within the scheme proposed 
by Grammaticos and Voros [20,21], result from an angular averaging in the 
momentum space. Studying the anisotropy in this momentum space might 
come out at large angular velocities as another source, and may be the main 
source, of discrepancy between the naive spin-independent liquid drop model 
approach and a fully self-consistent microscopic approach of pieces of nuclear 
matter experiencing fast collective rotations. 
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Fig. 1. Variation with the angular momentum / of three density parameters (po, 
a, 7) defined in the text, for the neutron distributions corresponding to four nuclei 
( 90 Zr, 150 Sm, 208 Pb, 240 Pu). The parameters a and po are expressed in fm and fm~ 3 
respectively. 




Fig. 2. Same as Fig. 1 for the proton distributions. 
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Fig. 3. Variation with the angular momentum / of the neutron {R n ) and proton 
(R p ) sharp density radius parameters (expressed in fm) corresponding to four nuclei 
( 90 Zr, 150 Sm, 208 Pb, 240 Pu). 
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Fig. 4. Root mean square energy deviation AE nuc (in MeV) resulting from the fits 
of the strong interaction part of the nuclear energy at various values of the angular 
momentum /, as a function of the number of nuclei A m i n defined in the text. 
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Fig. 5. Fitted values of the usual volume (a v ) and surface (a s ) energy coefficients 
expressed in MeV, as functions of the number of nuclei Np defined in the text for 
various values of the angular momentum /. 
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Fig. 7. Coefficients c v and c s denning the angular momentum dependence of the 
usual volume (a v ) and surface (a s ) energy coefficients as functions of the number of 
nuclei Np defined in the text. 
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Fig. 8. Root mean square radii of the proton distributions of four nuclei ( 90 Zr, 
150 Sm, 208 Pb, 240 Pu) expressed in fm as functions of the angular momentum /. 
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Fig. 9. Direct ("dir") and exchange ("exc") Coulomb energies of four nuclei ( 90 Zr, 
150 Sm, 208 Pb, 240 Pu) as functions of the angular momentum /. The l.h.s. (r.h.s. 
resp.) scale corresponds to the direct (exchange resp.) energies expressed in MeV. 




Fig. 10. Dependence of the results of the 2-parameters fit of the direct part of the 
liquid drop Coulomb energy in terms of the angular momentum /. Both parameters 
c and 7, defined in the text, are expressed in MeV. 
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Fig. 11. same as Fig. 10 for the 5-parameters fit. 
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Fig. 12. Dependence of the fitted value of the parameter c', expressed in MeV, of the 
liquid drop exchange part of the Coulomb energy defined in the text, as a function 
of the angular momentum /. 
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Fig. 13. Usual fissility parameter x of the U nucleus as a function of the angular 
momentum /. 
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